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The angular momentum operator in the Dirac equation 

Victor M. Villalbafl 



Abstract 

The Dirac equation in spherically symmetric fields is separated in two different tetrad 
frames. One is the standard cartesian (fixed) frame and the second one is the diagonal 
(rotating) frame. After separating variables in the Dirac equation in spherical coordinates, 
and solving the corresponding eingenvalues equations associated with the angular operators, 
we obtain that the spinor solution in the rotating frame can be expressed in terms of Jacobi 
polynomials, and it is related to the standard spherical harmonics, which are the basis 
solution of the angular momentum in the Cartesian tetrad, by a similarity transformation 
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1 Introduction 



The Dirac equation is a system of four coupled partial differential equation which describes 
the relativistic electron and other spin 1/2 particles. Despite the remarkable effort made 
during the last decades in order to find exact solutions for the relativistic electron in the 
presence of external fields, the amount of solvable problems is relatively scarce, being the 
Coulomb problem 0, the plane wave 0, and some electromagnetic configurations || the 
most representative examples. The problem of finding exact solutions of the Dirac equation 
is closely related to the possibility of accomplishing a complete separation of variables of 
the Dirac equation. 

During the last decades there have been different approaches to tackling the problem 
of separating variables in the Dirac equation in the presence of external fields; we have the 
Stackel spaces method developed by Bagrov and co-workers (|| based on the experienced 
acquired in separating variables in the Klein Gordon equation) , the theory of R-separability 
by Miller ||, |JJ , and matrix modification of the Stackel proposed by Cook || Recently an 
algebraic method of separation of variables [§, [HJhas been successfully applied in different 



gravitational and electromagnetic configurations. The main idea of the algebraic method 
of separation is to try to reduce the original Dirac equation, which in the presence of 
electromagnetic fields takes the form 

{D}V = {f(<9 M - T M - iAfj) + m} * = o (1.1) 
to a new equivalent equation according the following scheme 

{D}$> = ^{DjTT- 1 ^ = {Kx + K 2) <$> = (1.2) 

with 

^ = T$ (1.3) 
where Ki and K 2 are first-order commuting matrix differential operators 

'k Xl k 2 ] =0 (1.4) 



depending each one on different space-time variables, T is a constant nonsingular matrix, 
and JF that we have to find (if they exist) in order to fulfill the condition fll.4| ). De- 
pending on the problem to solve there are two different schemes of separation of variables 
(1) Consecutive separation of variables when one variable is separated from the remaining 
three variables, then one variable is separated from the other two variables, and finally, 
the remaining two variables are separated. (2) Pairwise separation of variables, when two 
space-time variables are separated from the another two ones, and the resulting operators 
are also separated. In the present article we focus our attention to the Dirac equation in the 
presence of a central field, like the Hydrogen Atom where a complete separation of variables 
is possible. In this case, using a pairwise scheme of separation, the angular dependence is 
completely separated from the other two variables. The article is structured as follows. 
In Sec. II we separate variables in the Dirac equation. In Sec III. we solve the angular 
equation in the local rotating frame in terms of Jacobi Polynomials In Sec. IV we establish 
the relation between the spinor solution in the diagonal (rotating) tetrad gauge and the 
fixed Cartesian one. 

2 Separation of variables 

In this section we proceed to separate variables in the Dirac equation ( [Lip , where 7 M are 
the curved gamma matrices satisfying the relation, 

{f\r} + = 2<r (2.i) 

and T M are the spin connections |7|] . 

T x = \g m [{dbl/dx x )a«-TZ x }s^ (2.2) 

with 

sr^rr-rr) (2.3) 

and the matrices , establish the connection between the Dirac matrices (7) on a 
curved spacetime and the Minkowski space (7) Dirac matrices as follows: 



% = b« la r = ojV (2.4) 

If we choose to work in the fixed Cartesian gauge, then spinor connections are zero and the 
7 matrices take the form 



7° = 7° = 7°, 7 1 



(7 1 cos ip + 7 2 sin ip) sin $ + 7 s cos •& 



7 



1 r 



7 



f 



(7 1 cos cp + 7 2 sin 99) cos $ — 7 s sin 1? 
-7 1 siny? + 7 2 cos<^) 



7_ 
r 



(2.5) 



73 



r sin r sin i? 

and the Dirac equation in the fixed tetrad frame (|2.5|) takes the form 



f{d t + iV(r) + fd r + ^ + —^- s d v + ml § = 

r r sin 1/ 



(2.6) 



In order to separate variables in the Dirac equation, we are going to work in the diagonal 
tetrad gauge where the gamma matrices 7^ take the form 



7° = 7°, f d = 7 1 , 7d = -7 2 , f d = ~^- q 

r r sin 17 



7 



(2.7) 



Since the curvilinear matrices 7 M matrices and 7^ satisfy the same anticommutation rela- 
tions, they are related by a similarity transformation, unique up to a factor. In the present 
case we choose this factor in order to eliminate the spin connections in the resulting Dirac 
equation. The transformation S can be written as || 



S 



1 



T7^ ex P(-75-7 1 7 2 )exp(-^7 3 7 1 )a = S a 



where a is the constant non singular matrix given by the expression 



(2.8) 



-(7V-7V + 7V + /) 



(2.9) 



which applied on the gamma's acts as follows 



37V 1 = 7 3 , 37V 1 = 7 1 , a7 3 a _1 = 7 2 , (2.10) 

the transformation S acts on the curvilinear 7 matrices, reducing them to the rotating 
diagonal gauge as follows 



S^YS = = 7^ (no summation) 



(2.11) 



then, the Dirac equation in spherical coordinates, with the radial potential V(r), in the 
local rotating frame reads 



L°d t + j% + +^d# + —^d v + m + z 7 V(r) 1 = 
[ r rsmv J 

where we have introduced the spinor related to \l/ by the expression 



(2.12) 



(2.13) 



and 7^ are the standard Dirac flat matrices. Applying the algebraic method of separation 
of variables ||, lQfl , it is possible to write eq. fl2.12|) as a sum of two first order linear 
differential operators K\ , K 2 satisfying the relation 



{k x + k 2 ) $ = o 



(2.14) 



KtQ = A$ = -K 2 <$> (2.15) 
then, if we separate the time and radial dependence from the angular one, we obtain 



7 2 ^ + 



smv 



7°7 1 $ 



(2.16) 



with 



k x ® = r U°d t + 7 1 <9 r + m + ryV] 7V $ = i«$ 



(2.17) 



$ = 7 ° 7 1 $ 



(2.18) 



where we have made the identification in = A. Notice that (|2.16|) is the angular momentum 
K obtained by Brill and Wheeler 0. Here it is necessary to remark that the operator K 2 
appearing in Q2.16|) is not singled valued and does not satisfy the properties of a "good" 



angular momentum operator. The true angular operator should be obtained from K 2 with 
the help of the transformation S ( J2.8D 

„3 



S 



l 2 d, + ^-d v 
sm v 



7 u 7 i 5- i 5$ = -ikS$ 



(2.19) 



since 



S 7 ^- 1 = 7^, we obtai 



am 



-i . 7 3 7°7 1 



7 7 7 «<? H + 7 7 7 SfyS H — - 

sm v sm 17 



■Sd in S 



-i 



(2.20) 



where S<£> = $ using the explicit form of the transformation S given by ( |2.8| ) we have 



Sd$S 1 = - cot i? H — (cos (p — / -f 1/ ~y 2 sin (p)^ 3 ^ 1 



(2.21) 



substituting ( ggj)) and (gjp into ( ^20]) 



(2.22) 



-0-2-1 <-> i 7 7' 7 r, ( 1 

7 7 7 ^ + . q 9 y + 7 U 
sin 17 



expression that can be written as follows 



ir = 7° (ctL)+1 



The operator K satisfies the relation 



K 2 = J 2 + \, J 2 = {L + \af 



and the eigenvalue k is related to j as follows 



(2.23) 



(2.24) 



(2.25) 



k 2 =j(j + l) + l = (j + \) 2 



(2.26) 



Now we proceed to decouple the equation ( |2.16| ) governing the angular dependence of the 
Dirac spinor In order to simplify the resulting equations, we choose to work with the 
auxiliary spinor $ related to $ as follows 



(2.27) 



then the equation (|2.16|) takes the form 



7% + -?—d v 
sin v 



7 ° 7 3 $ = _ m $ (2.28) 
In order to reduce the equation Q2.28Q to a system of ordinary differential equations, we 



choose to work in the following representation for the gamma matrices, |TT 



r 









l i J 




o 1 







Then, substituting ( |2.29| ) into ( |2.28| ) we obtain, 



a 



a 2 d$ — i -m 

sin v 



$i = — ik®i 



(2.29) 



(2.30) 



—a d$ + i— 



sin$ 



m 



with 



(2.31) 



$ 2 



(2.32) 



3 Solution of the angular equations 



In this section we are going to solve the systems of equations (|2.30|) and (|2.31| ). It is not 
difficult to see that the spinor $2 is proportional to <j 3 $i, this property allows us to consider 
only the system Q2.30I) . Using the standard Pauli matrices , we have that (|2.30|) reduces to 



fcsini^i — m^2 — sini?-^ = 

dv 



(3.1) 



k sin i9^2 — m£i + sin •& 



dii 

d-d 



(3.2) 



with 



6 

6 



(3.3) 



In order to solve the coupled system of equations (|3.1| )-( j3T2"D we make the following ansatz, 



6 = z(sin ^(cos ^) m+1 g(tf), 6 = (sin ^) m+1 (cos ^) m f($) (3.4) 
Substituting ( |3.4| ) into ( |3.1| ) and ( |3.2| ) we obtain, 

iJfeg(x)-(m+i)/(x) + (l-x)^^=0 (3.5) 

ikfix) + (m + J)g(:r) + (1 + x)^-^ = (3.6) 
2 ax 

where we have made the change of variable x = cos??, The normalizable solution of the 
system (|3.5|)-(P^D can be expressed in terms of Jacobi Polynomials Pj"'^ ]T2 



6 = C (sin -) m (cOS l)rn+lp(m-l/2,m+l/2) (3^ 

G = C (sin ^) m+1 (cOS l)mp(m+l/2,m-l/2) (gg) 

where Co is an arbitrary constant, which can be fixed in order to normalize the angular 
dependence of the spinor <3> according to the product 

2tt r dd = 6 kk , (3.9) 



o 



and n is given by the expression 



n = -m-- + k, n = 0,1,2.... (3.10) 



then, using the orthogonality relation for the Jacobi polynomials |~2] 



we obtain 



(1 - x) a (l + xfP^(x)P^(x)dx 

•1 

2 *+P+i Y(n + a + l)T{n + (3 + 1 



2n + a + p + 1 r(n + l)r(w + a + /3 + 1 



5mn (3-11) 



\/n\(n + 2m)! 

c = — '-= (3.12) 

(77 + 777,- 1/2)!^ V ' 



Now, we are going to solve the system of equation ( |2.17 ), governing the radial dependence 
of the spinor \1/ solution of the Dirac equation. This equation can be written in the form, 



^d t + 7 u <9 r + m 7 °7 3 - ir'V(r) + i- ) $ = 



K 



(3.13) 



Using the representation for the gamma matrices given by ( |2.29|) , we obtain the following 
system of equations 



d r + - J $1 - u yE - V(r) - m)$ 2 = 



(3.14) 



( -d r + -J $ 2 - a\E - V(r) + m)$i = 
therefore, the general form of the spinor $ is given by the expression 



■d i? 
$ = c (sin -) m (cos - 



m i(rrvp—Et) 



1 c(cosf)P,( m - 1 / 2 ' m+1 / 2 )A(r) N 
c(sinf)P^ +1 /2^-i/2)^( r ) 

(cosDP^- 1 / 2 ^ 1 / 2 )^^) 



V 



J 



(3.15) 



(3.16) 



where cq is given by ( p,12| ), c is a constant, and Air), B{r) satisfy the system of equations 



d r + - A(r) -(E- V(r) - m)B(r) = 
r / 



(3.17) 



-d r + - B(r) - (E — V(r) + m)A(r) = 



(3.18) 



obviously, the form A(r) and B(r) will depend on the potential V(r) . Among the solvable 
problems we have the Coulomb case where the radial dependence can be expressed in terms 
of confluent hypergeometric functions. [13] 

Regarding the eigenvalues m of the projection of the angular momentum operator —id v 
we have that since 



S g (<p + n) = -S g (<p) (3.19) 

and the spinor \l/ is single valued, then ^((p + 2tt) = and therefore m takes half 

integer values 



m = N + -, N = 0, ±1, ±2. 
2 



(3.20) 



4 Transformation of the spinor 



In this section we proceed to show that the angular dependence of the spinor solution of 
the Dirac equation ( |2.12| ), which is given by the expression ( |3.16| ), reduces to the standard 
spherical spinors when we apply the transformation S given by ( |2.8|) Since the spinor 
expressed in the Cartesian tetrad gauge, is related to the spinor solution \1/ in the local 
rotating frame by means of Q2.13D , we have 



5 , o7 7 3 $ = ^ 



(4.1) 



where we have used ( [2.181 ) an d (|2.10| ). Using the representation (|2.29|) and the spinor (|3.16| ), 
we obtain 



(cOs|)P ? ( m - 1 / 2 .-+l/2) fi ( r ) 

(sinf)P,( m+1 / 2 ' m ' 1 / 2 )5(r) 
-c(cosf)Pi m - 1 / 2 - m+1 / 2 )A(r) 
^ c(sinf)P ? ( m + 1 / 2 ' m ~ 1 / 2 )A(r) ) 



(4.2) 



e -i<p/2 














e iip/2 














e -i(p/2 














e i<p/2 



J 



cos f sin 2 

• «? ■& 
— sin ^ cos 2 













V 



o 
o 








cos | sin | 



sin g cos 2 



(4.3) 



substituting ( |4.3| ) into ( |4.2| ) and taking into account the relations 



X 



+ (l+x)P^ +1 \x) = 2P£ 



X) 



pw- 1 Hx) = p^ 1 0{x) + p<s!>{x) 



(4.4) 
(4.5) 



we find that the upper two components of the spinor \1/ 



take the form 



-/— — (sm - cos -J 1 



y/2r 

with the help of the relation 



(Sin | COS i)rn-l/2p(m-l/2,m-l/2)( CQS #y(m-l/2)<p B ^ 
I cos ^+l/2pKV^+l/2) (cos ^ e i (m+ l/2), jB(f) 



sin 



(4.6) 



1 /2Z + 1 



2 m l\ V 4vr 



(/ - m)\(l + m)!(l - x 2 ) m/2 P/™ m m) (x)e imv (4.7) 



and 



(1 - mV 

Pr m (x) = (-l) m 77-T3T^(^) 



(4.* 



(/ + m)\ 

we find that can be expressed in terms of spherical harmonics Yk,m{$, V 9 ) as follows 



= -i-(2k+ iy 1/2 
r 



(k + m+ l/2)y 2 Y™~ 1/2 ($, <p)B{ 



(4.9) 



^ (k-m + l/2)Y™ +1/2 {$, V )B{r) 
since the transformation So acting on the spinor \1/ has a block form, it is possible to write 
it as follows 



Sn 





r(sin , (9) 1 / 



then, if we set c=-l, we have 



cos I sm | 
- sm ~ cos | 



(4.10) 



and the lower components of the spinor \l/ read 



/ 



*~ = i-(2k + l)- 1/2 
r 



(k-m+ l/2) 1 / 2 Y™- 1/2 ($, <p)A(r) 



(4.11) 



(4.12) 



^ -(k + m + l/2)Y™ +1/2 ($, V )A(r) 
the results (|4.9|) and (|4.12 ) are in terms of the well known spherical spinors as in almost 
all the textbooks on relativistic quantum mechanics. The basis solution in terms of Jacobi 
Polynomials has been of utility in analyzing the Dirac spinor in the presence of a magnetic 
monopole, also it can be of help in discussing non-central interactions. 
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